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Employing the linear response Kubo formalism as implemented in a fully relativistic multiple- 
scattering Korringa-Kohn-Rostoker Green function method a systematic first-principles study based 
on density-functional theory (DFT) of the anomalous Hall conductivity (AHC) of the 3d-transition- 
metals Fe, Co and Ni is presented. To account for the temperature dependence of the AHC an 
alloy-analogy for a set of thermal lattice displacements acting as a scattering mechanism is used 
which is subsequently solved using the coherent potential approximation. Further, impurity scat- 
tering has been considered to elucidate the importance of an additional possible contribution to the 
AHC that might be present in experiment. The impact of correlations beyond the local spin-density 
approximation to the exchange-correlation functional in DFT is studied within the LSDA+f/ ap- 
proach. It is shown that both, the inclusion of correlations and thermal lattice vibrations, is needed 
to give a material-specific description of the AHC in transition-metals. 
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The simple experiment by Hall [l[ driving a current 
through a fcrromagnet and observing the anomalous Hall 
effect (AHE) as a transverse voltage has fruitfully spurred 
the development of experimental and theoretical meth- 
ods dealing with transport in solids. It now stands as a 
paradigm for understanding related transverse transport 
phenomena, like e. g. the spin Hall- (SHE), anomalous- 
and Spin-Ncrnst effect which have received intense in- 
terest in recent years. They all share a common origin, 
namely they are manifestly spin-orbit driven relativistic 
effects. 

The AHE has for decades eluded theoretical under- 
standing - it took more than 50 years until Karplus and 
Luttinger 0] put forward an insight which initiated mod- 
ern theories of the AHE. They identified the anomalous 
velocity as an interband matrix element of the current op- 
erator which is nowadays the foundation of semiclassical 
approaches which give a topological formulation of the 
AHE in terms of the Berry phase of Bloch bands in pure 
crystals [H, H| ■ The latter is used to define the so called 
intrinsic contribution to the AHE. Already early on Smit 
[H and Berger 0] discussed other extrinsic origins of the 
AHE, namely skew- and side-jump scattering. There are 
contributions to the AHE which fall in neither category, 
0] while it is now commonly established to separate the 
AHE into an intrinsic and a skew-scattering contribution 
and declare the difference to the total AHE as side-jump 
0. Experiments then rely on scaling mechanism to ex- 
tract these contributions from the raw data. 

Besides a wealth of model calculations (see the review 
0] and references therein) which arc tailored to iden- 
tify general trends but miss the material specific aspect 
a number of first-principles calculations building on a 
density-functional theory (DFT) framework employing 
the local spin density (LSDA) or generalized gradient 
(GGA) approximations have been undertaken recently 



to compute the anomalous Hall conductivity (AHC) in 
the transition-metals (TM) @-[lo|. Almost all of them 
rely on the Berry phase formulation for pure crystals and 
therefore are only able to deal with the intrinsic contri- 
bution. Boltzmann-transport theory based formulations 
have been used in the context of the SHE [ll| to compute 
the skew-scattering contributions in the dilute limit for 
alloys. Covering the whole concentration range of alloys 
and including all contributions to the AHE has recently 
been done flp . [T^ on the basis of a Kubo-Stfeda formu- 
lation 0, [ill . 

The role of correlations in the electronic structure of 
the 3d-TM has only very recently been addressed in the 
context of the AHE @,|loj]. Employing the LSDA turns 
out to give unfavorable agreement with experiment, and 
with the AHE being a property of the Fermi-surface 
[H [l(| it became clear that the LSDA does not sup- 
ply the proper bandstructure. This is demonstrated in 
particular for the case of Ni (see also table |T| where 
the LSDA/GGA strongly overestimates the magnitude of 
the AHC. Employing the LSDA/GGA+J7 remedies this 
problem, by moving down d-bands relative to the Fermi- 
energy (Ef), thereby making the X2 hole pocket present 
in LSDA/GGA disappear. 

A further important aspect of the AHE which is ad- 
dressed in experimental studies but rarely in theoreti- 
cal considerations is the temperature dependence of the 
AHE. For the pure 3<i-systems measurements of the AHE 
are typically done on commercially available specimen or 
thick layers grown on a substrate [17Hl9l | and the temper- 
ature is changed in order to vary the resistivity. The lat- 
ter makes the discussion of the temperature dependence 
very delicate when trying to disentangle different mecha- 
nisms and contributions to the AHE (inelastic scattering, 
scattering by phonons/magnons, etc.). It is advocated, 
however, as an empirical fact 0, [13] that inelastic scat- 
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tcring processes suppress the skew-scattering at higher 
temperatures with the intrinsic and side-jump (see how- 
ever remark above) contributions dominantly prevailing. 
This then again is used to experimentally analyze the 
AHE. Recently model calculations [20[ studied the role 
of inelastic scattering by phonons employing a Kubo- 
formalism and introducing a phenomenological scatter- 
ing rate 7 as the imaginary part of the self-energy. To 
our knowledge no first-principles approach has been used 
so far to deal with the temperature dependence of AHE 
in 3d-TM. 

In this letter we present a generally applicable formal- 
ism and results of a first-principles approach for calculat- 
ing the anomalous Hall conductivity of transition metals 
and their alloys. We show that the inclusion of both fi- 
nite temperature and correlation effects leads to a unified 
material specific description of these systems. 

As the AHE is inherently a relativistic phenomenon 
we choose to work within a fully relativistic approach 
employing the Kohn-Sham-Dirac equation as formulated 
in spin-polarized-DFT employing the Hamiltonian: 

H D = -ica ■ V + mc 2 f3 + V KS (r) + /3S ■ B xc (r) , (1) 

with Vks an d B xc being the spin-averaged and spin- 
dependent part of the one-particle potential, respectively, 
and the relativistic matrices ci, (3 and X having the usual 
meaning 21-23|. This has the important advantage [l5[ 
that disorder can be treated elegantly without making 
recourse to a Pauli approach which poses difficulties in 
calculating the vertex corrections. To determine longitu- 
dinal and transverse components of the conductivity ten- 
sor a natural starting point is the linear response Kubo 
framework which also can be used to derive the Berry 
phase related scmiclassical approach USEil. The Kubo 
approach has important advantages as compared to the 
latter. It allows straight-forwardly to include disorder, 
therefore not only being able to describe pure systems, 
but also alloys of the full concentration range including 
intrinsic and extrinsic contributions to the AHE 
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Further making use of an alloy analogy model (see be- 
low) finite temperatures can be accounted for. It also 
allows to include correlations beyond LSDA in the frame- 

26[. For cubic 
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work of LSDA+C/ or LSDA+DMFT 
and hexagonal systems with the magnetization pointing 
along the e 2 -direction, the AHE is given FJ, |l5| by the 
off-diagonal tensor element a yx = — <7 xy of: 

h 



Trace (j M (G+ - G~ 



)juG~ 



AniNil 



Trace ((G+ - G~ ) (r M j„ - r„^)) c ,(2) 



with the relativistic current operator j = — \e\ca and 
the electronic retarded and advanced Green functions G 
which in the framework of the presented KKR approach 



are given in a relativistic multiple scattering representa- 
tion [27j]. The angular brackets denote a configurational 
average which here is carried out using the coherent po- 
tential approximation (CPA) which allows to include ver- 
tex corrections (vc) , which are of utter importance for the 
quantitative determination of both the longitudinal and 
transversal conductivity in alloys. As has been argued 
and also shown 0, ED, [3] calculations omitting the vc's 
give the intrinsic AHC. Thereby, subtracting the latter 
from the AHC obtained from the value including the vc's 
the extrinsic part can be extracted. 

Several sources of electron scattering at finite temper- 
atures will determine the T-dependence of the AHE. We 
neglect the redistribution of states due to finite temper- 
ature in the electronic subsystem as well as electron- 
magnon interaction and consider as a dominant effect 
only thermal lattice vibrations. To include the latter as a 
source of electron scattering one could generalize Eq. @ 
to finite temperatures by including the electron-phonon 
self-energy S cl _ ph when calculating the Greens function 
G^ . This, however, is computationally very expensive. 
Therefore the consideration is restricted to elastic scat- 
tering processes by using a quasi-static representation of 
the thermal displacements of the atoms from their equi- 
librium positions as has already been used success fully 
by the authors in the theory of Gilbert-damping [28| . 
Treating each displaced atom as an alloy partner, we in- 
troduce an alloy-analogy model to average over a discrete 
set of displacements that is chosen to reproduce the ther- 
mal root mean square average displacement \J (m 2 )t for 
a given temperature T. This was chosen according to 
<« 2 >t = l^ffcl^i^ + \] with *(e D /D the De- 
bye function, h the Planck constant, k the Boltzmann 
constant and &d the Debye temperature [H|. Ignoring 
the zero temperature term 1/4 and assuming a frozen 
potential for the atoms, the situation can be dealt with 
in full analogy to the treatment of disordered alloys de- 
scribed above. To assess the presented approach the di- 
agonal resistivities (i. e. p xx ) were calculated for the 3d- 
TM making use of the Kubo-Greenwood expression for 
the symmetric part of the conductivity tensor (30l - [32l |. 
The results are shown as insets in Figs. [2] and [3] and com- 
pared with experimental data taken from the literature. 
As can be seen the agreement is rather good. Therefore 
we expect the aforementioned framework to be a reason- 
able approximation to properly describe electron-phonon 
scattering and the temperature dependence of the AHE. 

To study the impact of correlations (beyond LSDA) 
LSDA+J7 [§6| calculations have been performed keep- 
ing J off = 0.9eV fixed and scanning the [/-range up to 
typical values employed for the 3<i-TMs. In Fig. Q] the 
dependence of the intrinsic AHE at T = OK for bcc-Fe, 
hep-Co and fcc-Ni is shown. Whereas for Fe and Co only 
small variations of the AHC are observed, a pronounced 
{/-dependence for Ni is seen with the experimentally ex- 
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TABLE I. The intrinsic AHC (7 yx in (fi cm) 1 of the ferromagnetic transition metals Fe, Co and Ni from first-principles 
theoretical (present work compared to other) as well as experimental (Exp.) studies. The magnetization has been assumed to 
be oriented along the [001] direction. 
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FIG. 1. (Color online) The dependence of the AHC (T=0 K) 
bcc-Fe, hep-Co and fcc-Ni as a function of the {/-value in the 
LSDA+C/ calculation. 



tracted intrinsic value of -1100 S/cm recovered at a U- 
value of around 2.5eV (this value is also used in calcu- 
lation of the T-dependence below) . Analysis shows that 
this is due to a down shift of minority 3ti-bands w.r.t. 
Ep and a vanishing hole pocket at the X2 point as has 
already been recently discussed [l(| . In table U we show 
the calculated values for Fe (U = 1.8cV), Co (hep and 
fee, U = 2.3cV ) and Ni (U = 2.5eV) in comparison 
to other calculations as well as experiment. 

Early measurements of the AHE in Ni report a value 
of -646 S/cm at room temperature. [3^ Recent exper- 
imental work [i~8l | analyzed this in more detail claiming 
the AHE to consist of an intrinsic component of about 
-1100 S/cm and a sizable skew-scattering contribution at 
low temperatures, which both diminish at higher tem- 
peratures albeit with different rates. In Fig. [2] the cal- 
culated temperature dependence of the AHE in Ni using 
LSDA and LSDA+t/ (U = 2.5 eV, J e „ = 0.9 eV) as well 
as experimental results [l8| are shown. As could be ex- 
pected from the above the LSDA result strongly overesti- 
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FIG. 2. (Color online) The temperature dependence of the 
AHC of Ni. Shown are theoretical results obtained by LSDA, 
LSDA+(7 both including (vc) and excluding (nvc) vertex cor- 
rections and experimental data pl| . Triangle symbols denote 
the zero temperature values intrinsic AHE values for LSDA 
and LSDA+tT, respectively. The inset shows the longitudinal 
pxx component of the resistivity as compared to experiment 



mates the magnitude over the whole temperature range, 
whereas the LSDA+(7 fairly well reproduces the experi- 
mental result. This demonstrates that both correlations 
beyond LSDA as well as temperature induced thermal 
vibrations combined need to be taken into account. The 
vertex corrections due to the the lattice vibrations have 
little impact in the low T-regime and are negligible at 
higher temperatures such that, as seen in experiment, 
the intrinsic contribution survives. The deviation from 
the total AHC in the low T-range we attribute to possible 
impurities that might be present in the sample. In con- 
trast to Ni the temperature dependence in Fe (see Fig. [3]) 
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FIG. 3. (Color online) The temperature dependence of the 
AHC of Fe as calculated in LSDA and LSDA+[/. The in- 
sets show the longitudinal p xx component of the resistivity as 
compared to experiment [381 ]. 



FIG. 4. (Color online) The temperature dependence of the 
AHC of Nio.98Mgo.o2 (left) and Nio.gsFeo. 02 (right) as calcu- 
lated in LSDA+t/ compared to experiment [la ]. 



is found to be small. 

In the context of both the SHE and the AHE, 
39l . 4fj| it has already been shown that in the dilute/super 
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clean limit large skew scattering contributions can arise 
with the AHC scaling as er yx oc a xx . To demonstrate this 
and put it in the context of the recent experiment by 
Ye et al. [Lsj we performed calculation for Mg impuri- 
ties and Fe impurities in Ni. As can be seen in Fig [4] 
the calculated full AHC (including vc's) approaches the 
experimental curve for higher temperatures. However in 
the low temperature regime larger deviations are visible. 
Taking the difference between the calculation with vc's 
to those without vc's (intrinsic values), one also observes 
that the impurity induced extrinsic contribution for Mg 
shows the same sign as seen in experiment, i. e. it in- 
creases the absolute value, whereas for Fe impurities the 
opposite behavior is seen. This highlights again the fact 
that the skew-scattering component in an impurity spe- 
cific way determines the quantitative low temperature 
behavior of the AHC in clean 3c?-metals but also that 
the experimental determination of "clean" systems is ex- 
tremely challenging. 

In summary, we have used the linear response Kubo 
formalism as implemented in a fully relativistic multiple- 
scattering KKR Green function method to study system- 
atically from first-principles the anomalous Hall conduc- 
tivity (AHC) of the 3d-transition-metals Fe, Co and Ni. 
Going beyond the local spin-density approximation in 
DFT employing the LSDA+C7 and including finite tem- 
peratures by using a CPA-alloy analogy for the lattice 
displacements provided the necessary means to allow for 
a material specific description of the AHC. Further, the 
impact of dilute impurities has been analyzed. The pre- 
sented framework is now ready to be applied to the whole 
concentration range of correlated TM-alloys. Treating 
correlations beyond the static limit (LSDA+t/) of the 
LSDA+DMFT combined with a linear response trans- 



port formalism is a major issue for future work. 
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